Remote Sensing of Environment 140 17014: 60 77 



Contents lists available at SclcnccOirect 

Remote Sensing of Environment 

journal homepage: www elsevier corn/locatc/rse 



Evaluation of MODIS albedo product (MCD43A' over grassland, 
agriculture and forest surface types during dormant and 
snow-covered periods 



Cross Mark 


Zhuosen Wang Crystal B. Schaaf a . Alan H. Strahler b . Mark J. Chopping c , Miguel 0. Roman d . Yanmin Shuai cf . 
Curtis E. Woodcock b , David Y. Hollinger g , David R. Fitzjarrald b 

* Environmental. i'urtli and Ocvun Sciences Department, University n / Massachusetts Boston, Boston. MA USA 
*■ Center for Remote Sensing. Department ofCcography and Environment. Boston University. Boston. MA USA 
' Deitarlmenl of Barth and Environmental Studies. Montclair State University. Montclair. NJ. USA 
•’ Terrestrial Information Systems laboratory I Code 6 1 Of, NASA Coddard Space Flight Center. Creenbelt. MIX USA 
' Earth Ramtrres Tec hnology Inc,. Imirrl. MD. USA 
' NASA Coddard Space Flight Center, CreenbelL MD, USA 
8 USIM Forest Service Northern Research Station. Durham. NIL USA 
11 AtmospncnrScicncts Research Center, University ut Albany, SUNY, Albany, NV, USA 


ARTICLE INFO 


ABSTRACT 


Article history. 

Receiver! 25 October 2012 

Received in revised form 1 4 August 201 3 

Accepted 1 7 August 2013 

Available online 1 9 September 20 1 3 

Keywords: 

MODIS standard anti daily albedo pmducr 

Snow albedo 

forest 

Grassland 

Spatial representativeness 


This study assesses the Moderate-resolution Imaging Spcctroradiomeler (MODIS) BRDF/albedo 8 day slandard 
product and products from the daily Direct Broadcast BRDF/albedo algorithm, and shows that these products 
agree well with ground-based albedo measurements during the more difficult periods of vegetation dormancy 
and snow cover. Cropland, grassland, deciduous and coniferous forests are considered. Using an integrated 
validation slrategy, analyses of the representativeness of tire surface heterogeneity under both dormant and 
snow-covered situations are performed lo decide whether direct comparisons between ground measurements 
and 500-m satellite observations can be made or whether finer spatial resolution airborne or spjirbome data 
are required to scale the results at each location, landsat Enhanced Thematic Mapper Plus IETM f ) data are 
used to generate finer scale representations of albedo at each location to fully link ground data with satellite 
data. In general, results indicate the root mean square errors (RMSEs) are less than 0.030 over spatially represen- 
tative sites of agriculture/grassland during the dormant periods and less than 0.050 during the snow covered 
periods for MCD43A albedo products. For forest, the RMSEs are less than 0X120 during die dormant period and 
0.025 during the snow-coveted periods, However, a daily retrieval strategy is necessaty to capture ephemeral 
snow events or rapidly changing situations such as the spring snow inelL 

© 2013 Elsevier Inc All rights reserved. 


1. Introduction 

Surface albedo, defined as the ratio of the total (hemispheric) 
reflected solar radiation flux to the incidenL llux upon the surface, quan- 
tifies the radiation interaction between the atmosphere and the land 
surface. It plays a crucial role in land surface climate and biosphere 
models (Dickinson & Hanson. 1984; Dirmeyer & Shukla, 1994; 
Ijofgren, 1995). Carbon-only accounting approaches which ignore the 
albedo impacts of forests can significantly overestimate the climatic 
benefit of offsets (Thompson. Adams. 8? Johnson, 2009). As a key land 
physical parameter controlling the surface radiation energy budget, an 
absolute accuracy of 0.02 0.05 is required by climate models for global 
surface albedo (Dickinson. 1983. 1995; Dickinson et aL, 2008; 
Henderson-Sdlers & Wilson, 1983). 

• Corresponding author. TeL: I 1617287 5029. 

E-mail address: wangzhuoscnwginaitcom ( L Wang). 

0034-4257/S - sec front matter © 2013 Elsevier Inc All rights reserved 
http://dx.doi.org/ 1 0. 10 1 fi/j.rscZO 1 3.0S.025 


The surface reflectivity changes significantly with the appearance of 
snow. Snow also alters the exchange of moisture belween the surface 
and the atmosphere especially during the snowmelt period (Marshall. 
Roads. & Glalzmaier, 1994). The snow albedo contributes a strong 
positive feedback in surface modeling studies (Betts & Ball. 1997; 
Bonan, Chapin, & Thompson, 1995; Gardner & Sharp. 2010; Koltzow, 
2007; Li. Sun. Wang. & Liu, 2009; Molders. Luijting. & Sassen. 2008; 
Pedersen, Godtliehsen, & Roesch, 2008; Pedersen, Roeckner, Luthje, & 
Winther. 2009. Qu 8r Hall, 2006; Rutter et al.. 2009: Thomas & 
Rowntree, 1992; Viterbo 8i Betts, 1999; Wyser et al, 2007). The albedo 
at local solar noon for snow-covered forests is usually less than 03 in 
the shortwave, but it can reach 0.57 or higher for snow-covered grass- 
land and barren locations (Jin et al, 2002). 

The Moderate-resolution Imaging Spectroradiometer (MODIS) 
BRDF/albedo products (Schaaf, Wang, & Strahler. 2011a. Schaar el aL, 
2002. 2008) have been available since 2000 and provide high quality 
surface reflectance anisotropy retrievals over a variety of land surface 
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typos (jin ot jL 2003b; Liang el al„ 2002; Liu et .ll.. 2009; Roman ot al., 
2009, 2010], The Hidireri ino.il Reflectance Distribution Ftjnairm 
[BKDF) models and stupes (Gao. Scluaf. Strahler, Jin. & U 2003: Hiff 
Averin. Jian. SchaaT. A Amnion, 2008) art 1 also increasingly being used 
to provide information about the surface vegetation structure and the 
albedo quantities have been embraced by the global cltmale iruxleling 
and numerical assimilation communities (Dickinson et aL. 2008: Fang 
rt a I., 2007. Ijrwrence A Chase. 2007; Mairrelte, Barker, Cole, lacuna, 
A hncus. 2008: My lire. Kvalevag. A Scluaf. 200'j. LHeson et aL. 2003: 
Rnesdr, Scluaf, A Gao. 2004; Tun et al. 2004: Wei el aL 2001; Zliou 
et aL 2003). 

Allliougti reported on a 500 in grid, the MODLS BRDF/albedo 
prodtias an* usiully retrieved utilizing observations from a larger area 
depending an the view angles. Die footprint of observations increases 
with the increase in ihr scanning angfe (Tan et aL, 2006; Wolfe, Roy, A 
Vermote, 1998). Observauons are weighted by angular coverage before 
the retrieval bur it must still be acknowledged that the 500 m gridded 
product aaually represents information from a larger area. 

A number of studies luve validated the MODtS albrdu pmduii 
by a direct 'point-ro-ptxel* companson (Chen. Liang. Wang. Kim. A 
MartoiKhik, 2008 Liu et aL, 2009; Salomon. Schaaf, Strahler. Gao, A 
Jin. 2006; Wang et aL 2004) with the assumption that boih the tower 
and satellite albedo roughly sample the same spatial domain. However, 
ground measurements from tower usually only cover a limited area, 
which is much smaller than the MODIS spatial resolution. An important 
difficulty encountered whenever attempts are made to compare 
satellite-retrieved albedo values to tower- measured albedo daia is 
that the footprint of Ihe ground measurement is not always representa- 
tive of the greater satellite footprint Tftts is particularly true of rapidly 
changing surface conditions during senescence, green-op. and varying 
periods of snow cover. Assessments of spatial representativeness 
living Lindi.it data can pruvide an estimate of the general ability of 
the satellite retrievals to capture the tower measurements (Roman 
rt aL, 2009, 2010). However, conditions cm change so r.rpidly (luring 
dormancy and cplseincral snow events that such assessments must be 
reconsidered frequently for each situation. 

Initial investigations with "ponnt-to- pixel* comparisons at Ihe 
AmeriFlux Canadian University of California-lrvine (UCI) bum forest 
sites in Canada suggest rival rive MODIS winter values were usually 
lower than live lower values (Romdn et aL, 2009). However these loca- 
tions represent fickl sites with high spatial heterogeneity during winter 
periods. To consider the spatial scaling effect, finer spatial resolution 
I. inih.il linluncrd Thematic Mapper Plus (FTM r) data have been 
usually used to validate MODIS BRUF albedo products (Fang. Liang. 
Oven. Walthall, A Daughtry. 2004: Jin et aL 2003b liing et aL 2002; 
Susaki. Yasuoka. Kajiwara, Honda. & Hara. 2007). Barnsley et al 
(2000) aixl Luchl. Hyman. et al (2000' used Liixlsal data to anal>7r 
the spatial scaling effect on the albedo of a semi-desert environment 
prior to live launcii of MODIS. II ruuM lie ixiled however that must 
of these studies have assumed a lambcrtian surface to estimate 
landsat albedo although Jin el at. (2003b) calculated a Landval 
albedo Ivy applying the ratio of Ihr MODIS hemispherical albedo at 
local solar noon to the directional surface reflectance at the Landsat 
observing geometry. 

While the MODIS BRDF/albedo products have been evaluated during 
the growing season for a number of vegetated land covers with high ac- 
curacy (Disney. Lewis, Ihackralx Quaife. & Bamsiey. 2004; Jin et aL 
2003.1, 2003ft; Knisbelspiesse, Cairns, Schmid. Roman, A Schaaf. 2008; 
Lung et ai. 2002: Liu et aL 2009; Lynns. !m. A Randerson. 2008: 
Roescli et aL 2004; Roman et al.. 2009. 2010; Salomon el al.. 2006; 
Sanvain et al. 2008: Shuai. Scluaf. Shahlcr. Liu. A Jiaa 2008: Stone. 
Ancle i son, Sliettle. Aixlrews. A l.oukailmie, 2008; Susaki et aL 2007; 
Wang. Hartagc. Zeng. Dickinson. A Schaaf. 2005; Wang et al. 2004;. 
the assessment of the results during domunt and snow-covered sea- 
sons (Roman el al.. 2009. 2010: Schaaf. Liu. Gao, A Strahler. 2011b: 
WangetaL2012) bos only Just begun. 


The spatial patterns of albedo will change seasonally if live surface is 
comprised nf dilferent land types. This is especially tnie during dormanl 
periods (when vegetation is not photosyntheucatly active and leaves 
are either brown and/or very few). Jin rt aL (2002) analy/rd tlx- cried 
of snow over different land covers and showed tlut snow caused high 
heterogeneity in llie surface albedo making validation nxire diflicull. 
Therefore, differences between the tower measurements and MODIS al- 
bedo should be ex|tected and any evaluation of the MODIS BRDF/albedo 
products over dormant seasons, both snow-covered and snow-free, 
needs Id pay |tarticulai attention to spatial scaling effects. 

For agricultural and temperate grassland areas, the ground is often 
covered by dense uniform vegetation during the growing season. 
Thus, the surface is relative homogenous and ground measurements 
match welt with MODIS products (Chen el Jl, 2008: Jin et aL 2003b; 
liang el jl„ 2002; liu rt aL 2009; Roman ct al.. 2009, 2010; Susaki 
et al. 2007 j. However, these surfaces are often considerably mote liei- 
cTogenerxis during dnemant or partly snow -covered periods. Jin el aL 
(2003b; and Chen et aL (2008) showed that there are reiatrvely large 
differences Itel ween ground .iIIkhIos aixl MOWS alltrdn pnxluctv donng 
these seasons. 

Modeling and monitoring lire albedo of forests during the snow 
covered periods can be difficult. A snow-covered forest in winter is 
made lip an uppei layer of dark leaves and branches and a lower layet 
of bright snow on the forest floor. Because of the upper tbliage in 
dense evergreen forests, the reriected photons from the snow surface 
have difficulties escaping the forest canopy especially undei high zenith 
viewing angles, moreover, a large fraction ol the uixferstory may lie 
shaded in winter scenes, owing to the higher solar zenith angles 
However, the reflectance from deciduous forests, where snow informa- 
tion is nxire visible through a canopy of bare branches, is relatively high, 
rise density of the canopy can also affect whether a pixel is designated 
as snow-covered or snow- free. 

This study aims to assess tlx* accuracy of MODIS albedo products by 
cnrnpiroixi with grouixl measurements after establishing Ihe spatial 
representativeness of vegetated surfaces during dormant and snow- 
covered situations Sectxin 2 describes both tire ground measurements 
and MODIS alhrxlo products and Section 3 outlines the assessment stra- 
tegy. In Section 4, we compare tire ground measurements to MODIS 
albedos and lands.it albedos where appropriate and discuss the overall 
accuracy of MODtS products during rite domunt and snow-covered 
seasons. 


2. Datasets 

2.1. Ground measurements 

The primary validation sacs used for the MODtS Albedo product have 
historically been the Surface Radiation Budget Network (SURFRAD) sites. 
(Augustine. Drfuisi. A long. 200(1) which are maintained in the United 
States by NOAA as part of the Baseline Suffice Radiation Network 
(BSRN) (Ohmura rt ai. 1998). Hu* World Climate Research Programme 
(WCRP; Radiative Fluxes Working Group initiated rise BSRN to support 
the research projccls or ihe WCRP aixl other sriemifx' programs. These 
sewn SURFRAD sites (Fort I’eck. Nri. Sioux Foils. SD, IVnn State, PA. 
Boodville, ll. Table Muunl.un, CO. Goodwin Creek, MS and Dr-setl Rock. 
NV] (Table 1 ) ore used again in tlus study. In addition to the seven 
SURFRAD sites, the Boulder, CD site is also another BSRN site dose to 
the Table Mountain site that utilizes pyranometers mounted on a 
300 ni tower. These etgltf sites |xuvide rise higlsesl quality inlenalilxattsl 
albedo measurements. 

Tluee additiixi.il sites of New Frigiond lorests from llie AmeriFTux 
network are also used in this study. AmcriRux was established by the 
Department ol Riser gy (DOF.) in 1996 to provide continuous otaervu 
tions of crosysicm level exchanges of OOz. water, and energy, including 
surface albedo (Law et al. 2002 Running et aL. 1999). 
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TaMr I 

SUWKAIl RSRN. and Ameriflux ground sites. 


Matinri name 

Latitudc/longitude 

tower hcight/footprinl (m] 

Slate 

land cover 

MOOCS tile 

Nrtwock 

Konkin 

4UU5/ 

I0S01 

300/378825 

Colorado 

Grassland 

1KUVIM 

1ISKN 

fort Pock 

4031/ 

105.10 

10/126.28 

Montana 

Grassland 

HI 1 VD4 

SURfRAD 

Goodwin Crock 

34.25/ 

R9R7 

10/12628 

Mississippi 

Crassljnd 

Hinvra 

NUKIKAD 

Sioux Falls 

4173/ 

9662 

10/12028 

South Dakota 

Grassland 

HI 1 VW 

SURfRAD 

Table Mkjunutr. 

ICl'l7ix28 

Colorado 

Grassland 

II09VIM 

SUklKAU 


4QI3/ 
ins. 74 
Desert Rock 

3662/ 

1 1 11(13 

1002028 

Nevada 

Desert, sparse grass 

HCHVtB 

SURfRAD 

BonctviUc 

4IIIIV 

mi/ 

UVUtx/H 

lllinon 

Agnculhitc 

imviM 

MIKIKAI) 

Penn Sun- 

4072 / 

7733 

100202a 

PcniKytvanii 

AsrictiHuie 

HIZVtM 

SURfRAD 

Harvard EMS 

«3W 

72.172 

30060 

Massachusetts 

Mixed forest 

nt 2 vw 

AmmHux 

lltm land larch 

4S216' 

68.709 

30066 

Maine 

ncnx'ntf-% ruvdIHeaf fort-sl 

H13V04 

AnirriRux 

llontand West 

4520EV 

18.747 

30966 

Maine 

Ewiixeen needleleaf forest 

IlUVtM 

AmeriHux 


Total downward and upward radiation (028-3.0 cm) is measured by 
Eppley pyranonve-ters mounted on the 10 m towers at the SURFRAD 
sites while the Boulder site utilizes a 300 m tower. Normal Incidence 
Pyraiwmeters and shaded pyranometers are used to measure the direct 
normal and diffuse shortwave fluxes (Augustine et aL. 2000). The 
instruments are ventilated and heated during the winter so that there is 
a great deal of confidence in the winter measurements as wvll as the 
growing season ones. Kipp and Zonen albedometers in the shortwave 
(0.3-221 cm) arc used to measure albedo at AmeriFlux sites These 
albedometers ate mounted on 30 m towers and not routinely (seated 
and ventilated, then-fore the data are not considered as reliable during 
snow precipitation periods. Ground albedo is calculated by the ratio of up- 
welling radiation and downwelling radiation during die local mid-day 
time. The footprints of 10 m. 30 m. and 300 m towers are estimated to 
be about 126 m. 366 m and 3788 m in diameter respectively. 

/ n 2H can HFOV 1 

where /is the circular footprint of ground tower measurements. H |m| is 
the tower height, and HFOV |degrees| is its half of field of view. HPOV is 
81° (Michalsky, Harrison. & Berkheiser, 1995). 

the Bondville SURFRAD site is located southwest of Champaign. 
Illinois, ui an agricultural region ( Fig. 1 ). This site represents a mixture 
of crops and drainage ditches which are maintained with a variety of 
harvesting and fallowing practices. It is quite heterogeneous at the 
MODES scale (Liu el al„ 2009; Salomon et al. 2006). The Sioux Falls 
SURFRAD site is located on die grounds of the Earth Resources Observa- 
tion and Science (EROS) Data Center, South Dakota, and is covered by 
grass. The Goodwin Creek site, west of Oxford. Mississippi, ts located 
on rural pasture land surrounded by deciduous trees. The Penn State 
University (PSU) SURFRAD site ts located on the grounds of PSU's 
agricultural research farm about 6 miles southwest of State College, 
Pennsylvania and is in a broad Appalachian valley between theTussey 
and Bald Eagle Ridges. The Fort Peck SURFRAD station is located on 
the Port Peck tribes Reservation in Montana. This site is dommated by 
native grasses. The Table Mountain SURFRAD station is located in 
Colorado and is a mix of exposed rocks, sparse grasses, desert shrubs 
and small cactus. The Boulder BSKN station, also in Colorado, is very 
close in the Table Mountain station and is covered by grass. Tlx- Desert 
Rock SURFRAD suuon is covered with very sparse vegetation and lies to 
the ixirt hwest of Las Vegas, Nevada. The AnieriFlux site at I lowland For- 
est ts located in central Maine. USA. about 35 miles north of Bangor. 
Tlx- natural stands in this boreal-northern hardwood transitional forest 
consist of hemlock-sproce-fir. and hemlock -hardwood mixtures. The 
evergreen needle-leaf canopy height at Howland West Tower is about 
20 m. The albedo data used here arc from 2007 to 2009, and the 
alliedometer was Ixsiled during the winter of 2009. The Howland 
larch Tower is surrounded by deciduous needle-leaf trees (Lvch. lorix 
kinchin] (Hollinger. Ollinger, Richardson. Meyers, Ik Dail, 2010). The 


albedo data were collected in 2008 for the Larch Tower. The AmeriFlux 
Harvard Environmental Measurement Station (EMS) tower at Harvard 
Forest lies in the central Massachusetts town of IVtcrsham. The domi- 
nant species indude red oak. red maple, black birch, white pine, and 
hemlock. 

22. Satellite albedo 

22 l MCD43A standard HKI)l : /albedn product 

Tlie standard MODIS BRDF/albedo product (MCD43A) (lucht, 
Schaaf. & Strahler, 2000; Schaaf et al„ 2002] provides the weighting pa- 
rameters associated with the RossThick-LiSparse Reciprocal (RTLSS; 
BRDF model that describes the reflectance anisotropy of each pixel al 
a 50O-m gndded resolution. Both a sufficient number of observations 
and good angular sampling are needed to achieve a full inversion 
retrieval and estimate a high quality BRDF. Acknowledged as a tradeoff 
between the temporal stability of the surface reflectance and tlw avail- 
ability of sufficient angular samples, a 16-day period of cloud-free, 
atmospherically-corrected surface reflectances from both Terra ( MOD) 
and Aqua (MYD) is used to derive MCD43A BRDF/albedo (Gao. Schaaf. 
Strahler. ft Lucht. 2001; Roy. Lewis. Schaaf, Devadsga. ft Boschelli, 
2006; Wanner ct al., 1997). Wilh an 8 day system of overlapping pro- 
cessing. more phenological variability can be accurately characterized. 
A backup algorithm (also called a magnitude inversion) is employed if 
a lugh quality full inversion retrieval cannot be accomplished due to 
poor angular sampling or insufficieni input observations (Schaar el al.. 
2002). This a priori knowledge method often performs quite well 
under normal situations (Jin et al. 2003a, 2003b; Liu et al, 2009; 
Salomon et at 2006) but should be considered a poor quality result 
and is assigned a poor quality flag. The MCD43A BRDF/albedo standard 
product only retrieves a snow albedo quantity when the majority of 
observations during a 16 day period are snow -covered. Snow-covered 
and non-snow observations are currently always processed separately. 

The intrinsic black-sky albedo (BSA) at local solar noun and tlx? 
white-sky albedo ( WSA) is generated by integrating the BRDF calculat- 
ed from the three retrieved parameters (/»> / geo and f val ). Blue sky 
albedo, which considers both the diffuse and direct incident radiation 
for a specific time and atmosplieric suite, can be calculated as follows 
(Lewis & Barnsley, 1994): 


Muc-dcy 1 SKYI. , ^c»sl> 1 SKYL j lUark-iky I 


where SKYlf i) b the proportion of diffuse irradiation at a certain solar 
zenith angle ,. The SKYI is measured by the shaded pyranometcrat thc 
SURFRAD and Howland Forest (West and Larch) sites. Haivjrd FMS 
Forest does not have ground SKYI. data and the MODIS aerosol product 
(MOODS) is used to calculate the proportion of diffuse irradiation. 
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SHF RAD. BSRN, and 
AmeriFlux Sites from 
Google map 



Fig.1 . The locations of the eleven SUFRAD, BSRN and AmeriFlux sites Boulder, C.oodwin Creek, Fort Peck, Sioux Falls, Table Mountain. Desert Rock, Bondvlllc, Penn State, Harvard EMS. 
Howland West and Howland larch sites from Coogle map 


In addition to the spectral quantities, three broadband albedos of 
VIS (0.3-0.7am), NIR (0.7-5.0OH) and shortwave (0.3-5.0im; are 
calculated from the 7 spectral bands via narrow to broadband 
conversion coefficients (Liang, Strahler, & Walthall, 1999; Stroeve 
etal., 2005;. 

shortwave snowiree 03973 , 02382 2 0 3489 3 0 2655 , 

0:1604 5 0 0138 6 0 0682 7 0 0036 

shortwave snow 0:1574 , 02789 2 03829 3 0:1131 5 

00694 7 0 0093 : 4 

22.2. MCI Ml A direct broadcast ( DR) albedo data 

The MODIS DB BRDF/albedo product (Shuai, 2010) which is utilized 
for regional applications (e.g., forest, agriculture and disturbance moni- 
toring) is operated in a daily rolling mode to provide more frequent 


surface Nadir BRDF-Adjusted Reflectances (NBAR) and albedos than 
the current standard product. Two versions of the DB mode are utilized 
in this study. First the 16-day daily mode uses 16 days worth of reflec- 
tances as input but weights the more recent clear observations 
(with the highest observation coverage) more heavily. A retrieval is 
attempted each day based on the proceeding 16 days' worth of data. 
Full inversions are performed if there are sufficient high quality and 
well sampled angular observations, where the quality of the angular 
sampling is determined by reference to the Weights of Determination 
(Lucht & lewis, 2000). Otherwise a magnitude inversion will be 
processed with the a priori backup database being updated wiLh the 
most recent full inversion results. Second, the 1-day daily mode, a 
magnitude inversion is performed each day using the full inveisions 
from the 16-day daily mode as the a priori information for each 
succeeding day. While this mode emphasizes the single day observation 
the most heavily, it is also more sensitive to any cloud contamination or 
residual aerosols affecting that single look. 



